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Abstract. We study the synchrotron emission of neu- 
trino pairs by relativistic, degenerate electrons in strong 
magnetic fields. Particular attention is paid to the case 
in which the dominant contribution comes from one or 
several lowest cyclotron harmonics. Calculations are per- 
formed using the exact quantum formalism and the qua- 
siclassical approach. Simple analytic fits to the neutrino 
synchrotron emissivity are obtained in the domain of mag- 
netized, degenerate and relativistic electron gas provided 
the electrons populate either many Landau levels or the 
ground level alone. The significance of the neutrino syn- 
chrotron energy losses in the interiors of cooling neutron 
stars is discussed. 



1. Introduction 

In this article, we study the synchrotron emission of neu- 
trino pairs by relativistic, degenerate electrons in strong 
magnetic fields, 



(1) 



This emission can be an important source of neutrino en- 
ergy losses from the crusts of magnetized cooling neutron 
stars (NSs) as well as (in a modified form) from the su- 
perfluid NS cores. 

The process has been considered by several authors 
starting from the pioneering article by Landstreet (1967). 
The adequate quasiclassical formalism was developed by 
Kaminker, Levenfish & Yakovlev (1991, hereafter KLY) 
who presented also critical analysis of preceding work. 
KLY carried out detailed analytic and numerical analysis 
of the case in which the electrons populated many Lan- 
dau levels and the main contribution into the synchrotron 
emission came from high cyclotron harmonics. They also 
considered briefly the opposite case of the superstrong 
magnetic field which suppresses greatly the contribution 



from all harmonics. In that case, only the first harmonics 
actually survives (although exponentially damped). 

Exact quantum formalism of the neutrino synchrotron 
emission was developed by Kaminker et al. (1992a). 
Kaminker & Yakovlev (1993) considered the synchrotron 
process in the nondegenerate electron gas. Recently 
Kaminker et al. (1997) have analyzed how the synchrotron 
emission is modified in the superconducting NS cores 
where the initially uniform magnetic field splits into flux- 
oids. In a recent article, Vidaurre et al. (1995) have recon- 
sidered the results of KLY claiming them to be inaccurate. 

In the present article, we continue to study the neu- 
trino synchrotron radiation from the degenerate magne- 
tized electron gas. The emphasis is made upon the case in 
which the electron transitions associated with one or sev- 
eral lowest cyclotron harmonics are most important. This 
case has not been studied attentively earlier. We calcu- 
late the neutrino emissivity numerically using the exact 
quantum formalism (Sect. 2, Appendix A). We also use 
the quasiclassical approach (Sect. 3, Appendix B), com- 
bine our new results with those by KLY and propose an 
analytic fit that describes accurately all the neutrino syn- 
chrotron radiation regimes in a strongly degenerate, rela- 
tivistic electron gas in which the electrons occupy many 
Landau levels. Another limiting case in which the electrons 
populate the ground level alone is considered in Appendix 
C. In Sect. 3 we present also critical discussion of the re- 
sults by Vidaurre et al. (1995). In Sect. 4 we show the 
importance of the neutrino synchrotron emission in the 
NS interiors. 

2. Quantum formalism 

We will mainly use the units in which m e = c = h = 
&b = lj where fee is the Boltzmann constant. We will 
return to ordinary physical units whenever necessary. The 
general expression for the neutrino synchrotron energy loss 
rate (emissivity, ergs s _1 cm -3 ) from an electron gas of 
any degeneracy and relativity, immersed in a quantizing 
magnetic field, was obtained by Kaminker et al. (1992a): 
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x Au;f(l-n (2) 

Here, Gp = 1.436 x 10 -49 ergs cm 3 is the Fermi weak- 
coupling constant and b = B / B c is the dimcnsionless mag- 
netic field (B c = m 2 e c 3 /(he) w 4.414 x 10 13 G). Further- 
more, n and p z are, respectively, the Landau level number 
and the momentum along the magnetic field for an elec- 
tron before a neutrino-pair emission; the energy of this 
electron is e = yjl + 2nb + p 2 . The primed quantities n' 
and p' z refer to an electron after the emission; its energy is 
e' = ^l + 2n'b + p' 2 ; f = /(e) = {cxp[(e - p)/T] + l}- 1 
is Fermi-Dirac distribution of the initial-state electrons, 
/' = f(e') is the same for the final-state electrons; p is 
the electron chemical potential and T is the temperature. 
The energy and momentum carried away by a neutrino- 
pair arc denoted as to = e — e' and q, respectively. The 
z component of q is q z = p z — p' z , while the component 
of q across the magnetic field is denoted by q±. The sum- 
mation and integration in (2) is over all allowed electron 
transitions. The integration has to be done over the kine- 
matically allowed domain q 2 + q\ < uo 2 . The differential 
transition rate A (summed over initial and final electron 
spin states) is given by Eq. (17) of Kaminker et al. (1992a). 
Taking into account that some terms are odd functions of 
p z and vanish after the integration, we can rewrite A in a 
simple form, 

A = { ^ ql q ^ {pl +P ' 1 + 2 ) + ^ ( * $) 
- (u 2 -ql-ql) 2 ^+(u 2 -ql-ql)<$>) 

~ ^[(2^-2^ -£)(*-*) 

+ 3 (c 2 - ql q\) *] . (3) 

Here, p± = \/2nb and p' ± = \]2n'b are the transverse 
momenta of the initial-state and final-state electrons, re- 
spectively; 

* = ^-i.W + ^.-iW, 

^ = ^'-l,n-lW + ^',nWi (4) 

u = q 2 J(2b), F n , n {u) = (-l) n '- n F nn ,(u) = 

u (™-™')/2 e -"/2( n M/ n !)l/2 jL »-«'( u ) ; and is an as _ 

sociated Laguerre polynomial (see, e.g., Kaminker & 
Yakovlev 1981). Furthermore, in Eq. (3) we introduce 
G\ = E v (Cv + C a) ^d C 2 _ = ^(Cy ~ C 2 A ), where C v 
and Ca are the vector and axial vector weak interaction 
constants, respectively, and summation is over all neutrino 
flavors. For the emission of the electron neutrinos (charged 
+ neutral currents), one has C'v = 2sin 2 #w + 0.5 and 
Ca = 0.5, while for the emission of the muonic or tauonic 



neutrinos (neutral currents only), C' v = 2sin 2 #w ~ 0.5 
and C' A = —0.5; #w is the Weinberg angle. Adopting 
sin 2 6» w ^ 0.23 we obtain C\ w 1.675 and C 2 w 0.175. A 
comparison of the neutrino synchrotron emission by elec- 
trons with the familiar electromagnetic synchrotron emis- 
sion has been done by KLY. Although electromagnetic 
radiation is much more intense it does not emerge from 
deep neutron star layers (neutron star interior is opaque 
to photons) while neutrino emission escapes freely from 
stellar interior, producing an efficient internal cooling. 

We have composed a computer code which calculates 
Qsyn from Eqs. (2) and (3) for arbitrary plasma param- 
eters pY e , T and B, where p is the mass density and Y e 
is the number of electrons per baryon. Technical details 
are presented in Appendix A, the results are illustrated 
in Sect. 4, and the case of very strong magnetic field, in 
which the bulk of electrons occupy the ground Landau 
level, is considered in Appendix C. 

3. Quasiclassical treatment 

In this section, we develop quasiclassical description of the 
neutrino synchrotron emission from a degenerate, ultrarel- 
ativistic electron gas in a strong magnetic field B — 10 11 - 
10 14 G. Therefore we assume that n ^> m e c 2 and T -C Tp, 
where Tp = [fi — ra e c 2 )jkp, w /u/fce is the degener- 
acy temperature. We mainly analyze the case in which 
the electrons populate many Landau levels. This is so if 
the magnetic field is nonquantizing or weakly quantizing 
(e.g., Kaminker & Yakovlev, 1994), i.e., if fj, » \f\ + 2b. 
The latter condition is realized at sufficiently high den- 
sities p > p B , where p B « 2.07 x l0 6 b 3 / 2 /Y e g cm~ 3 = 
2.23 x 10 5 B% 2 /Y e g cm" 3 , and B 13 = B/(10 13 G). At 
these densities, p, is nearly the same as without magnetic 
field, p w epp w c?i(37r 2 n e ) 1 / 3 , where pp is the field-free 
Fermi momentum. The formulated conditions are typical 
for the neutron star crusts. 

In the case of many populated Landau levels, we can 
replace the summation over n by the integration over p± in 
Eq. (2). The remaining summation over n' can be conve- 
niently replaced by the summation over discrete cyclotron 
harmonics s — n — n'=l, 2,3,... Furthermore, an initial- 
state electron can be described by its quasiclassical mo- 
mentum p and pitch-angle 9 (p z = pcos8, p± — psinO). 
Since we consider strongly degenerate electrons we can set 
e = p and p = pp in all smooth functions under the inte- 
grals. Then the integration over p is done analytically. In 
this way we transform the rigorous quantum formalism of 
Sect. 2 to the quasiclassical approximation used by KLY 
(see their Eqs. (3), (4) and (8)). 

According to KLY the quasiclassical neutrino syn- 
chrotron emission is different in three temperature do- 
mains A, B, and C separated by two typical temperatures 
T P and T B - 

T P = = It b x 3 « 2.02 x l0 9 B 13 x 2 K, 

2fc B 2 
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Tn = 



Tvuj* b 
k B 



1.34 x 10 9 Bi 3 (l +X 2 )- 1 ' 2 K. 



(5) 



Here, x = p-p/(m e c) is the relativistic parameter (x ^> 1, 
in our case), and u> B = eBcj p, is the electron gyrofre- 
quency at the Fermi surface. Notice, that while apply- 
ing our results to the neutron star crusts, one can use 
a simplified expression x w 1.009(^6 ^e) 1 ^ 3 , where p e = 
p/(10 6 g cm" 3 ). 
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Fig. 1. Density-temperature domains of different neutrino 
synchrotron regimes. Tp, is temperature (5) below which (do- 
main C) the emission goes via the ground cyclotron harmon- 
ics s = 1 (for B = 10 12 , 10 13 and 10 14 G, logB shown in 
parenthesis). Tp is temperature (5) below which (in domain 
B) the Pauli principle restricts the number of cyclotron har- 
monics (see text); it is shown only for B = 10 12 G because Tp 
becomes higher than the electron degeneracy temperature T F 
for B = 10 13 and 10 14 G. Domain A (T P < T < T F ) is real- 
ized only for B = 10 12 G; Tf is independent of B, for displayed 
parameters. 



Figure 1 demonstrates the main parameter domains 
A, B and C for the ground-state (cold catalyzed) matter 
of NS crusts. Thermal effects on the nuclear composition 
are neglected which is justified as long as T < 5 x 10 9 
K (e.g., Hacnsel et al. 1996). It is assumed that nuclei 
of one species are available at any fixed density (pres- 
sure). Then the increase of density (pressure) is accom- 
panied by jumps of nuclear composition (e.g., Haensel et 
al. 1996). The ground-state matter in the outer NS crust, 
at densities below the neutron drip density (4 x 10 11 g 



cm -3 ), is described using the results by Haensel & Pichon 
(1994) based on new laboratory measurements of masses 
of neutron-rich nuclei. At higher densities, in the inner 
NS crust, we use the results of Negele & Vautherin (1973) 
derived on the basis of a modified Hartree-Fock method. 
Small discontinuities of the curves in Fig. 1 are due to the 
jumps of the nuclear composition. Notice, that the prop- 
erties of the neutrino synchrotron emission vary rather 
smoothly in the transition regions from domain A to B 
and from B to C. 

The high-temperature domain A is defined as Tp <C 
T <C Tf; it is realized for not too high densities and 
magnetic fields where Tp <C Tf- In Fig. 1, this domain 
exists only for B = 10 12 G at p < 10 11 g cm~ 3 . In do- 
main A, the degenerate electrons emit neutrinos through 
many cyclotron harmonics; typical harmonics is s ~ x 3 . 
Corresponding neutrino energies u> ~ lo b x 3 -C T are not 
restricted by the Pauli principle. The quasiclassical ap- 
proach of KLY yields 

2 G F k B Tm 2 e uj B x s 



Q. 



A 

syn 



1897T 5 c^Tl 

3.09 x 10 15 Bf 3 T 9 x s 



(25C 2 - 21C±) 



ergs cm 3 s , 



(6) 



where Tg = T/(10 9 K). Here and in what follows, numeri- 
cal factors in the practical expressions are slightly different 
from those presented by KLY because now we use more 
accurate value of the Fermi constant (see Eq. (2)). 

The moderate-temperature domain B is defined as 
Tb < T <C Tp and T <C Tf- It covers wide temperature 
and density ranges (Fig. 1) most important for applica- 
tions. In this domain, neutrinos are again emitted through 
many cyclotron harmonics s ~ k^T /Tilu* b 3> 1, but their 
spectrum is restricted by the Pauli principle, and typical 
neutrino energies are to ~ ksT. As shown by KLY, in this 
case the neutrino emissivity is remarkably independent of 
the electron number density: 



2C(5) G 2 F m 2 e Lo 2 B 
9tt 5 c 5 h 8 
9.04 x 10 14 B 2 3 T| ergs cm" 3 s _1 , 



G\ (k B T) 5 



(7) 



where C(5) w 1.037 is the value of the Riemann zeta func- 
tion. 

The third, low-temperature domain C corresponds to 
temperatures T < Tp at which the main contribution into 
the neutrino synchrotron emission comes from a few lower 
cyclotron harmonics s=l, 2, . . . If T <C Tp, even the first 
harmonics s = 1 appears to be exponentially suppressed 
as discussed by KLY. A more detailed analysis will be 
given below. 

The emissivity in the combined domain A + B, 
including a smooth transition from A to B at T ~ Tp, 
was calculated accurately in KLY. The results (Eqs. (13), 
(15) and (18) in KLY), valid at Tp < T < T F , can be 
conveniently rewritten as 



(8) 
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Sab 



27£* 



7T 2 2 9 C(5) 
Tp 3 



r 



f + (0-^f_(0 



(9) 



where the analytic fits to the functions -F±(£) read 
(1 + ciyi) 2 



F-(0 = D 2 



(l + am + bnjf)^' 

1 + c 2 y 2 + d 2 yl + e 2 y| 
(1 + a 2 y 2 + &22/I) 5 



(10) 



13/2 



In this case y 1>2 = [(1 + ai i2 $ 2/3 ) 2/3 • 
1CT 4 , bt = 7.405 x 10~ 8 , ci = 3.675 x 10~ 4 , 
1CT 3 , 6 2 = 1.536 x 10~ 5 , c 2 = 1.436 x 10~ 2 
10" 5 , e 2 = 7.647 x 10~ 8 , D x = 44.01, D 2 
3172, a 2 = 172.2. 

Now let us study the neutrino emissivity in the com- 
bined domain B + C. Using the quasiclassical expressions 
of KLY in the ultrarelativistic limit (m e — > 0) we obtain 



ai = 2.036 x 

d 2 = 1.024 x 
36.97, ai = 



)BC _ ,qB c 
' syn ^fcsyn^ 



BC, 



(11) 



where 



S 



1 



BC 



with 



2 6 C(5) z 2 T 7 



J sin 3 6» d6» J q ± dq A 



J dq z (u 2 - q 2 z - - $)- 



* - $ = 2 



1 J^) + J^) 



Here, J s {y) is a Bessel function of argument y = 
P±<l±c/{eB), and J^(y) = dJ s (y)/dy. The neutrino-pair 



energy can be expressed as w k g z cos 



The 



tegration should be done over the kinematically allowed 



domain: 



> 



sin 6(n z - q z f 



Kj_ = SU) B 



I sin 9, 



k z = suj* b cos 9/ sin 2 9. 

One can easily see that Sbc depends on the only pa- 
rameter z. The domain B corresponds to x~ 3 <C z -C 1 
while the domain C corresponds to z > 1. Equation (12), 
in which we set m e — > 0, does not reproduce the high- 
temperature domain ^4. However, we have already de- 
scribed the transition from B to .4 by Eqs. (8-10). 

At z <C 1, the neutrino emissivity, determined by 5bc 
in Eq. (12), comes from many cyclotron harmonics (see 
KLY, for details). A Bessel function J s (y) and its deriva- 
tive can be replaced by McDonald functions. The main 
contribution to the integrals (12) comes from a narrow 
vicinity — q±\ < cj^s 1 / 3 , and \q z — k z \ < lu b s 2 ^ 3 of the 
saddle point q_\_ = kj_ and q z — n z , in which the neutrino- 
pair energy is nearly constant, uj fa sui B / sin 2 9. Adopt- 
ing these approximations and replacing the sum over s 



by the integral we reproduce evidently the result of KLY, 
Sbc — 1 • However if wc take into account small corrections 
in the expression of a Bessel function through McDonald 
functions (Sokolov & Ternov 1974), and weak variation of 
u near the saddle-point, we obtain a more accurate asymp- 
tote Sbc = 1 — 0.4535 z 2 / 3 . The derivation is outlined in 
Appendix B. 

In the opposite limit of z > 1, one can keep the con- 
tribution from the first harmonics s = 1, and replace 
(* - $) -> 1, (e"/ T - l)- 1 -> e" w / T in Eq. (12) as 
described in KLY. At any kinematically allowed q z and 
qx, the neutrino-pair emission is suppressed by a small 
factor ~ e _w / T . As shown in KLY, the most efficient neu- 
trino synchrotron radiation occurs from a small vicinity of 



the allowed region, where q z w — 
mum and the neutrino energy to ; 



■ cost 



has 



} / (1 + cos 9) is most 



strongly reduced by the quantum recoil effect. This region 
corresponds to the backward electron scattering. The in- 
tegration over q± is then done analytically, and we are left 
with a two-fold integration over q z and 9. In the limit of 
very high z, it gives 



Sbc 



2C(5) 



exp 



(-1) 



1 



28 



(14) 



The convergence of this asymptote is very slow, and we 
present the second-order correction term 28/z which im- 
proves the convergence considerably. For instance, at z = 
60 the two-term asymptote gives an error of about 6.5%, 

(12) 

while the one-term asymptote gives an error of about 36%. 

Notice, that the emissivity given by Eq. (20) in KLY 
in the limit of z ^ 1 is 4 times smaller than the cor- 
rect emissivity presented here (due to a simple omission 
(13) in evaluating Q^ yn made by KLY). Thus Eqs. (20) and 
(21) in KLY are inaccurate at z ^s> 1. 

In addition to analyzing the asymptotes, we have cal- 
culated Sbc numerically from Eq. (12) in the quasiclassi- 
cal approximation at intermediate z. The results are fitted 
by the analytic expression 



'BC 



exp(-z/2)D 1 (z)/D 2 (z), 



(15) 



where L>i(z) = 1 + 0.4228z + 0.1014z 2 + 0.006240 z 3 , 
D 2 (z) = 1 + 0.4535 z 2 / 3 + 0.03008 z - 0.05043 z 2 + 
0.004314 z 3 . The fit reproduces both the low-z and the 
high-z asymptotes. The rms fit error at z < 70 is about 
1.6%, and the maximum error is 5% at z sa 18. 

Now, we can easily combine Eqs. (8) and (|ll]) and ob- 
tain a general fit expression for the neutrino synchrotron 
emissivity which is valid everywhere in domains A, B, C 
(T <C Tp, p 3> Pb), where the electrons are degenerate, 
relativistic and populate many Landau levels: 



Q. 



ABC 



Qsyn Sab Sbc- 



(16) 



Here, Q^, n is given by Eq. (7), while Sab and Sbc are 
defined by Eqs. (9), (10) and @. 
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B 13 

Fig. 2. Neutrino synchrotron emissivity vs B from a plasma 
with pY e = 10 11 g cm -3 and T = 10 9 K. Solid curve shows our 
quasiclassical result, dash-curve — calculation of Vidaurre et 
al. (1995), dots — analytic Eq. (26) of these authors. 



The neutrino synchrotron emission in domains B + C 
(in our notations), has been reconsidered recently by Vi- 
daurre et al. (1995). Their results are compared with ours 
in Fig. 2 which shows the neutrino synchrotron emissiv- 
ity as a function of B for the same plasma parameters 
(pY e = 10 11 g cm~ 3 , T = 10 9 K) as in Fig. 4 by Vi- 
daurre et al. (1995). For the parameters chosen (due to 
very high p), domain A is not realized, while the tran- 
sition from domain B to C occurs at fairly high B (for 
instance, z — 1 corresponds to B « 3.5 x 10 14 G). Our 
quasiclassical calculation, analytic fit (|l^) and quantum 
calculation are so close that yield the same (solid) line. 
Numerical calculation of Vidaurre et al. (1995) is shown 
by dashes, and their analytical approximation (claimed to 
be accurate at intermediate z) by dots. In our notations, 
the latter approximation reads S^} — 0.0073 z 5 e~ z . We 
see that such an approximation reproduces neither low-z 
nor high-z asymptote, and is, in fact, inaccurate at in- 
termediate z. It also disagrees with numerical curve of 
Vidaurre et al. (1995) and with our curve. The numerical 
results by Vidaurre et al. (1995) are considerably differ- 
ent from ours except at B sa (6 — 8) x 10 15 G. At lower 
and higher B the disagreement becomes substantial. In 
domain B, Vidaurre et al. (1995) present another analytic 
expression that differs from Eq. (7) by a numerical fac- 
tor 9/(27r). The difference comes from two inaccuracies 
made by Vidaurre et al. (1995). First, while deriving the 
emissivity, they use the asymptote of the McDonald func- 
tion at large argument (their Eq. (19)) instead of exact 
expression for the McDonald function (as was done by 



KLY) . One can easily verify that this inaccuracy yields an 
extra factor 3/7T. Secondly, Vidaurre et al. (1995) calcu- 
lated inaccurately an integral over q z (their Eq. (22)) as 
if the integrand were independent of q z . This yields the 
second extra factor 3/2. Thus we can conclude that the 
results by Vidaurre et al. (1995) are rather inaccurate in a 
wide parameter range. We have performed extensive com- 
parison of the results obtained with our quantum code 
and with the quasiclassical approach. We have found very 
good agreement (within 3-5 %) in all the cases in which 
the quasiclassical approach can be used (see above). This 
statement is illustrated in Figs. || and |J. Figure [| shows 
the quantum and quasiclassical synchrotron emissivities 
from a plasma with pY e = 10 8 g cm -3 and T = 10 9 K as 
a function of B. Figure || displays the ratio of the quan- 
tum to quasiclassical emissivities in more detail. The qua- 
siclassical approach is valid as long as p ^ ps (as long 
as electrons populate many Landau levels) which corre- 
sponds to B < 6 x 10 14 G for our particular parameters. 
For these magnetic fields, the quantum and quasiclassi- 
cal results are seen to coincide quite well. Low-amplitude 
oscillations of the curves in Fig. |] reflect oscillations of 
the synchrotron emissivity as calculated with the quan- 
tum code. The oscillations are produced by depopulation 
of higher Landau levels with increasing B. They represent 
quantum effect associated with square-root singularities 
of the Landau states. As seen from Fig. ^ the oscillations 
are smeared out with increasing T by thermal broadening 
of the square-root singularities ( cf. Yakovlev & Kaminker 
1994). In a higher field B > 6 x 10 14 G, the electrons 
populate the ground Landau level alone and the electron 
chemical potential is reduced by the magnetic field. Very 
high B, not shown in Fig. |§|, remove the electron degener- 
acy. Note that at high temperatures, one should take into 
account the synchrotron neutrino emission by positrons 
(Kaminker & Yakovlev 1994). All these conditions are de- 
scribed by the quantum code while the quasiclassical ap- 
proach is no longer valid. In Appendix C, we obtain simple 
asymptotic expressions for Q syn in the case in which the 
electrons populate the ground Landau level alone. Cor- 
responding curves are shown by dashed curves in Fig. Bl 
and they are seen to reproduce the exact quantum curves 
quite accurately. 

4. Discussion and results 

Figures ||-|| display density dependence of the neutrino 
synchrotron emissivity Q syn calculated from Eq. ( fL6| ) for 
the magnetic fields B = 10 12 , 10 13 , and 10 14 G at four 
temperatures T = 3 x 10 9 , T = 10 9 , 3 x 10 8 , and 10 8 K, 
respectively. We adopt the ground-state model of matter 
in the NS crust (see Sect. 3). Various neutrino-emission 
regimes can be understood by comparison with Fig. 1. 

The high-density (horizontal) parts of the synchrotron 
curves correspond to domain B (Eq. (7)), where Q syn is 
density independent. The low-density bends are associated 



6 



V.G. Bezchastnov et al.: Neutrino synchrotron emission 



at 

g 

a 



O ° 



t = i n 9 k 


J_ JA 1U J- V 




— _ T - Sxld K 

SN 


\ 

^\ \ 
\ \ \ 




v \ 
vx \ 

VS. \ 




v s \ 
v \ \ 




\ \ \ 




\ \ \ 




\ \ V 




\ \ \ 


pY e — 10 8 g cm 3 









100 



B 13 



Fig. 3. Neutrino synchrotron emissivity versus B from a 
plasma with pY e = 10 8 gem" 3 at T = 3 x 10 8 , 5 x 10 8 
and 10 9 K. Solid curves are calculated with the quantum code, 
long-dash-curves are obtained using the quasiclassical calcula- 
tion valid at B < 6 x 10 14 G, dash-curves show the quantum 
asymptotic expression (Appendix C) for the case in which the 
ground Landau level is populated alone, B > 6 x 10 G. 
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Fig. 5. Density dependence of neutrino emissivities from 
ground-state matter of the NS crust due to various mechanisms 
at T = 3 x 10 9 K. Curves 'syn' '(12)', '(13)', and '(14)' refer to 
the synchrotron mechanism at B = 10 12 , 10 13 and 10 14 G, re- 
spectively. Curve 'pairs' corresponds to neutrino emission due 
to annihilation of electron-positron pairs; it is almost indepen- 
dent of B at given T. Other curves are for B = 0: 'brems' — 
the total electron-nucleus bremsstrahlung; 'plasma' — plas- 
mon decay; 'photo' — photoneutrino process (see text). 




0.8 1 1 1 1 1 1 1 1 1 

15 20 25 30 35 40 45 

B13 

Fig. 4. Ratios of the quantum to quasiclassical neutrino syn- 
chrotron emissivities presented in Fig. H. 




(14) - 
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log p [gem ] 

Fig. 6. Same as in Fig. ^ but at T = 10 9 K. Pair annihilation 
depends noticeably on B. 
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T=3 x 10 K 




9 10 11 12 13 

log p [g cm ] 

Fig. 7. Same as in Figs. | and | but at T = 3 x 10 s K. Pair 
annihilation and photoneutrino processes become negligible. 
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log p [g cm ] 

Fig. 8. Same as in Figs. | - ^ but at T = 10 s K. Plasmon 
decay becomes negligible. 



with transitions either into domain A (where Q syn cx p 8 ^ 3 
according to Eq. (6)) or into domain C (where Q syn de- 
creases exponentially due to cyclotron harmonics suppres- 
sion, Eq. (E3)). Domain A is realized only for B = 10 12 G 
and T > 10 9 K in Figs. || and ||. The low-density bends 
°f Qsyn in domain C are much steeper than those in do- 
main A. These bends are more pronounced at highest 
B = 10 14 G, at which domain C extends to higher T and 
p (Figs. | and |). 



For comparison, we also plot the emissivities produced 
by other neutrino generation mechanisms: the electron- 
positron pair annihilation into neutrino pairs, electron- 
nucleus bremsstrahlung, plasmon decay and photon decay. 
The pair annihilation in a magnetized plasma has been 
considered by Kaminker et al. (1992a, b), and Kaminker 
& Yakovlev (1994). For the parameters of study, the emis- 
sivity appears to be weakly dependent on the magnetic 
field. At B < 10 13 G it is very close to the zero-field emis- 
sivity (Itoh et al. 1989, 1996). As seen from Figs. || and 
§, the pair-annihilation emissivity differs slightly from the 
zero-field one only in a not too hot plasma (T < 10 9 K) 
at B > 10 14 G. The neutrino bremsstrahlung curves are 
plotted neglecting the influence of the magnetic field. The 
effect of the field on the bremsstrahlung has not been stud- 
ied so far but it is expected to be weak, for the param- 
eters in Figs. &-§[ We use the results of Haensel et al. 
(1996) to describe the neutrino pair bremsstrahlung due 
to Coulomb scattering of electrons by atomic nuclei in the 
liquid phase of matter. In the solid phase, similar process 
is known to consist of two parts: the phonon and static 
lattice contributions. We use the results by Yakovlev & 
Kaminker (1996) to evaluate the phonon contribution. As 
for the static lattice contribution, we employ the most re- 
cent theory by Pethick & Thorsson (1996) and perform nu- 
merical calculation from Eqs. (28) and (29) of their paper 
(adopting the Debye- Waller factor and the nuclear form- 
factor which were used by Yakovlev & Kaminker 1996). 
Numerous jumps of the bremsstrahlung curves in Figs. 
H are associated either with jump-like changes of nuclear 
composition of cold-catalyzed matter or with solid-liquid 
phase transitions (see Haensel et al. 1996 for details). The 
neutrino emissivities from other processes are determined 
by the electron and positron number densities which are 
nearly continuous function of the density. Therefore, all 
other curves are smooth. The neutrino generation due to 
plasmon and photon decays in a magnetic field has not 
been considered in the literature, and we present the field- 
free results of Itoh et al. (1989, 1996), for illustration. 

In the case of zero magnetic field, the bremsstrahlung 
process dominates completely in most dense layers of the 
NS crust at not too high temperatures T < 10 9 K. Plas- 
mon decay, photon decay, and pair annihilation are signif- 
icant at high temperatures, T > 10 9 K, but their emissivi- 
ties become negligible very soon as temperature decreases. 

The synchrotron emissivity is, to some extent, similar 
to the bremsstrahlung, for it persists over the wide tem- 
perature and density ranges. In the presence of the strong 
magnetic field B > 10 13 G, the synchrotron emission is 
seen to be important and even dominant for any T in 
Figs. |]-||. In a hot plasma (Fig. ||), the synchrotron emis- 
sion is significant at comparatively low densities, p = 10 s - 
10 9 g cm' 3 . With decreasing T the neutrino synchrotron 
emission becomes more important at higher densities. At 
T = 10 8 K, only the bremsstrahlung and synchrotron 
emissions actually survive (Fig. §h; if B = 10 14 G, the 



8 



V.G. Bezchastnov et al.: Neutrino synchrotron emission 



synchrotron emission dominates over the bremsstrahlung 
one in a wide density range, p > 10 9 g cm -3 . 

5. Conclusions 

We have considered the synchrotron neutrino-pair emis- 
sion of electrons from a dense magnetized plasma. We de- 
veloped a computer code (Sect. 2) which calculates the 
synchrotron emission in the presence of quantizing mag- 
netic fields for a wide range of conditions at which the elec- 
trons can be weakly as well as strongly degenerate and/or 
relativistic. We have also calculated the synchrotron emis- 
sivity in the quasiclassical approximation (Sect. 3) for 
strongly degenerate, relativistic electrons which populate 
many Landau levels. We have paid special attention to the 
case in which the main contribution into the synchrotron 
neutrino emission comes from the fundamental or several 
low cyclotron harmonics. This case has not been analyzed 
properly earlier by KLY and Vidaurre et al. (1995). We 
have obtained a simple analytic expression ( |l6| ) which fits 
accurately our quasiclassical results in wide ranges of den- 
sities, temperatures, and magnetic fields. In Sect. 4 we 
have demonstrated that the synchrotron neutrino emis- 
sivity gives considerable or even dominant contribution 
into the neutrino emissivity at B > 10 13 G in moderate- 
density and/or high-density layers of the NS crusts, de- 
pending on temperature and magnetic field. Let us no- 
tice, that the neutrino emissivities plotted in Figs. 5-8 
are appropriate for cooling neutron stars rather than for 
newly born or merging neutron stars. Since the neutrino 
radiation from a cooling neutron star is too weak to be 
detected with modern neutrino observatories we do not 
calculate the spectrum of emitted neutrinos. 

Our results indicate that the neutrino synchrotron 
emission should be taken into account in the cooling theo- 
ries of magnetized neutron stars. It can be important dur- 
ing initial cooling phase (t— 10-1000 yrs after the neutron 
star birth). At this stage, the thermal relaxation of inter- 
nal stellar layers is not achieved (e.g., Nomoto & Tsuruta 
1987, Lattimer et al. 1994) and local emissivity from differ- 
ent crustal layers can affect this relaxation and observable 
surface temperature. 

The synchrotron process in the crust can also be sig- 
nificant at the late neutrino cooling stages (t ~ 10 5 yrs, 
T ~ 10 s K), in the presence of strong magnetic fields 
B > 10 14 G. The synchrotron emission can be the domi- 
nant neutrino production mechanism in the neutron star 
crust, while the neutrino luminosity from the stellar core 
may be not too high, at this cooling stage, and quite com- 
parable with the luminosity from the crust. 

It is worthwhile to mention that the synchrotron emis- 
sion can be important also in the superfluid neutron star 
cores (Kaminker et al. 1997). A strong superfluidity of 
neutrons and protons suppresses greatly (e.g., Yakovlev 
& Levenfish 1995) the traditional neutrino production 
mechanisms such as Urea processes or nucleon— nucleon 



bremsstrahlung. The superfluidity (superconductivity) of 
protons splits an initial, locally uniform magnetic field 
of the neutron star core into fluxoids — thin magnetic 
threads of quantized magnetic flux. This process modifies 
the neutrino synchrotron process (Kaminker et al. 1997) 
amplifying it just after the superconductivity onset and 
making it significant since the traditional neutrino gener- 
ation mechanisms are suppressed. It has been shown that 
this modified neutrino synchrotron process can dominate 
over other mechanisms if the initially uniform magnetic 
field in the NSs core is B > 10 13 G and if T < 5 x 10 8 K. 
These fields and temperatures are quite consistent with 
the results of the present article: similar conditions are 
necessary for the synchrotron emission to dominate over 
other neutrino production mechanisms in the NS crust by 
the end of the neutrino cooling stage. If so, the total neu- 
trino luminosity (from the stellar crust and core) can be 
governed by internal stellar magnetic fields which can af- 
fect the neutron star cooling. We plan to study this cooling 
in a future article. 
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Appendix A 

For performing exact quantum calculation of the neutrino 
synchrotron emissivity from Eq. (2) we replace the in- 
tegrations over p z and q z by the integrations over ini- 
tial and final electron energies e and e' , respectively. In 
this way we can accurately integrate within the intervals 
|e — fi\ < T and \e' — fi\ < T, where the Fermi-Dirac dis- 
tributions are rapidly varying. It is also convenient to set 
q± = (uj 2 — q^j (1 — rj) and replace the integration over q± 
by the integration over rj. Then 

G 2 b » roo d£ 
s=l J V £ ~ £ n> 

D = [ drj {C\R+ - C 2 _R-) , 
Jo 

R + = {1 + r, [l+pl+p'l - (cc 2 - ql) 77] } 
- [1 + (P±+P'l)r)]$, 

R- = (1 + rj)V- (1 - 2??)$. (Al) 
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Here, n' = n — s (s enumerates cyclotron harmonics); 
e n — vl-|- 2nb and £ n < = \/l + 2n'b correspond to the 
excitation (de-excitation) thresholds of the initial and fi- 
nal electron states, respectively. The neutrino pair energy 
and longitudinal momenta are determined from conserva- 
tion laws: p z = yje 1 — e 2 and p' z — ±y/ ' e a — e 2 n , . The 
sign of p' z and the domain of integration over e' are speci- 
fied by the minimum longitudinal momentum of the final- 
state electron for given n, e and s: 



P'l < P'z < Pi, 



Pl,2 



2(sTPz) 



(A2) 



If p[ > 0, one has p' z = y/ e' 2 — e 2 , and e\ < e' < e' 2 , 
where 



4' + {STPZ? 

2{eTPz) 



(A3) 



If p'i < 0, then there are two integration domains. The 
first one is e n i < e' < e[, with p' z = — \Je' 2 — e 2 , . The 
second domain is e n > < s' < e' 2 , with p' z = \J s 12 — et,. 

Appendix B 

Consider the z — > asymptote of the function Sbc( z ) 
given by Eq. (12). The integrand of (12) contains the fac- 
tor R = w 2 /(exp(w/T) — 1). Since the neutrino-pair energy 
w = q z cosO + su>* B varies near the saddle point, we have 



R(ui) = R{lu q ) + R' Aoj + -R" Aw 2 , 



(Bl) 



where lu q — k z + suj* b 
Aw = cos8(q z — k z ) , 

2 + e"(-4 + Av- 



is the value of w in the saddle point, 



R" 



'(2 



(e*> - 1)3 



(B2) 



v = cjq/T, and n z = suj b cos 8/ sin 2 8 is the z-coordinate 
of the saddle point. The term in ( |Bl| ), which is linear in 
Aw, does not contribute to (12) due to integration over 

1z- 

In the low-z limit (i.e., for T ^> Tb), the Bessel func- 
tions J s (x) and J' 8 (x) in jl^ ) can be expressed through 
McDonald functions K-L/z(rj) and ^2/3(^7) of the argument 
j] = (s/3)e 3 / 2 , where e = 1 - (x/s) 2 = 1 - {q ± /K±_) 2 : 
k±_ = lu b s/ sin 8 is the coordinate of the saddle point 
transverse to the magnetic field. These expressions can 
be written as (Sokolov & Ternov 1974) 



J s (x) = 
J'six) = 



ttV3 
e 

e 
5 



#1/3 



2/3 



2K 



2/3 



Jq (K1/3 - 6vK 2 / 3 



1 

3?7 



3r7 K 1/3 



(B3) 



KLY took into account the main terms of these asymp- 
totes. Here, we include small corrections (proportional to 
the factor e < s~ 2 / 3 <C 1 in square brackets). This yields 



2c 2 
3tt 2 



/3 



K; 



2^1/3-^2/3 



2/3 
1 



4^ 2 2 /3 



6ry 



(B4) 



Let us substitute (glj) and (B4) into (12). Then 

Sbc(z) = 1 + a lZ 2/3 + a 2 z 2/3 = 1 - 0.4535z 2/3 . (B5) 

In this case, a\ = LxM\Ni/ (10 Iq) = 0.2053 comes from 
variation of w, with 

d6 sin 5/3 8 cos 2 8 = 0.45890, 

dvR"v 10/3 = 3 2/3 x 76.533, 

poo 

/ dr 1 ^ 3 {K 2 1/3 + K 2 2/3 ) = 0.76706, 
Jo 

(B6) 

Furthermore, a 2 = L 2 M 2 N 2 /(5 Iq) = —0.6588 comes from 
the corrections to the McDonald functions, 







Li 


= [' 




Jo 


Mi 


= 3 2 / 3 




/'OC 


m 






Jo 


h 







L 2 
M 2 
N 2 





32/3 



d8 sin 5/3 6» = 1.6826, 
„io/3 



dv 

e*-l 



20.468, 



drj rf 



8/3 



6K 2 1/3 4 



4K, 



2/3 



-2.6082. 



(B7) 



Appendix C 



Consider the neutrino synchrotron radiation by relativis- 
tic degenerate electrons in a very strong magnetic field 
(6^1, vl + 26 > fx) in which the bulk of electrons pop- 
ulate the ground Landau level. The electron Fermi mo- 
mentum is then given by pp — 2a; 3 / (36), and the chemical 
potential is fi = y/l + p\ . The number of electrons on the 
excited Landau levels is exponentially small, and the ma- 
jor contribution to the neutrino emission comes from the 
electron transitions from the first excited to the ground 
Landau level. Equation (3) reduces to 



C 2 

A = — [(1 + 6 - u Q b)(2bu - 2bu + u) - u] e~ u , 
ee' 



(CI) 



where uq — (lu 2 — q 2 )/(2b). Inserting ( |Cl"| ) into (2) we can 
integrate over q± = \j2bu: 



h 2 C 2 

Vsyn — / 



3(2tt) e 



d Pz I dq z uf(l-f) 

>{u >0} 
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x — {(1 + b - buo)[2bu + 1 - 26 

EE 1 

+ (2b - 1 - uoje-™ ] + (1 + w )e" 1 



1}. (C2) 



Since the first Landau level is almost empty, we can set 
/ « exp(-(e-/i)/i), where t = T/(m e c 2 ) = T/(5.93 x 10 9 
K)). The neutrino emission is greatly suppressed by the 
combination of Fermi-Dirac distributions /(l — /'). These 
distributions determine a very narrow integration domain 
which contributes to Q sy n- 

Equation (C2) can be further simplified in the two JJ 
limiting cases. The first case is fi c < // < yl + 26, 
or pf/2 < 6 < b c , where fj, c = (1 + b)/y/l + 26 and 
6 C = Pf{pf + m)- The main contribution into Q sy n comes u\ 
from narrow vicinities of two equivalent saddle points 
Pz0 = ±[(pF + ^) 2 -26-l]/[2(p F + ^)], q z0 = =FV(PF + M). 
Each point corresponds to the most efficient electron tran- 
sition in which an initial-state electron descends to the 
ground Landau level just with the Fermi energy s' = fx, 
emitting a neutrino-pair with the energy ui = | q zQ | . In this 
case, the energy of the initial-state electron is e = 
One has uq <C 1 in the vicinities of the saddle points. 
Expanding e, s', and uq in these vicinities in powers of 
(pz —Pzo) and (q z — q z o), we obtain from Eq. (C2) 



Q 



G| C\b A 



(26 + 3)i 4 



syn 



3(27T) 4 p|.(pF + /i)-H' 3 5 



exp 



(C3) 



where H = I — (6/6 c ) and 5 = sin(7ri7). 

The second case corresponds to /i < /i c or 6 > b c . Now 
the most efficient electron transitions are those in which 
the initial-state electron is near the bottom of the first 
Landau level (e ps \/l + 26, ~ + 26). Accordingly, 
w = Vl + 26 - e' and u = (1 + 6 - eVl + 26)/6. The 



energy of the final-state electron e' 
from the maximum energy e' = (1 



varies mostly 



b) /VI + 26 (asso- 
ciated with the minimum allowable neutrino-pair energy 
a; = 6/Vl + 26 at uo = 0) to the minimum energy e' = fi 
allowed by the Pauli principle. One can put p z = in 
all smooth functions under the integral (C2), and replace 
integration over q z by integration over e' or over u, with 
(1 — /') = 1 m the integration domain. The integration is 
then taken analytically (for 6^1) yielding 



Qsyn — 



G|C246 5 a 



3(2tt) 5 (1 + 26) 3 / 4 
x F(ui) exp 



VT+26 
t 



where 



F{u x ) = ^--ui + 2-2e" 



— u\ e 



(C4) 



(C5) 



and U\ = 1 — (fi/ /j, c ). In the limit of [i — > 1 we have U\ — * 1 
and F = (4/3) - (3/e)). Then Eq. (C4) reproduces the 
asymptotic expression for the synchrotron neutrino emis- 
sivity of electrons from a non-degenerate electron-positron 



plasma (Eq. (35) of the paper by Kaminkcr & Yakovlev 
1993). In the opposite limit, /i — > /i c , we obtain u\ — * 
and F fa uf/6. 

Both asymptotes, ( |C3] ) and (C4), become invalid in a 
narrow vicinity |1 - (p/fi c )\ « 0.5|1 - (6/6 c )| < V*/& 1/4 
of the point p = fi c or 6 = b c . We propose to extend 
the asymptotes, somewhat arbitrarily, to the very point 
fi = fi c by replacing 




7T7, 



(C6) 



in Eqs. (C3) and Eq. (C4). These replacements do not af- 
fect significantly Qsyn outside the vicinity of /i = /i c but 
produce physically reasonable interpolation within this 
vicinity. By matching the modified asymptotes at fi = ji c 
and implying again > 1 we get 7 = 1.905t 1 / 2 6 _1/ ' 4 . In 
this way we obtain a complete set of equations to describe 
Q syn at /j, < VI + 26. 
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